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The distribution function of transmitted charge through a 
double-barrier junction is studied at zero temperature and at 
small applied voltage. Both a semiclassical model, in which 
the transport is described by jump rates, and a quantum 
mechanical model, which averages over resonant and non- 
resonant states, are used to determine the characteristic func- 
tion of the transmitted electrons. It is demonstrated that 
for large times the logarithm of the characteristic function is 
equal within the two approaches. The charge distribution is 
in between a Gaussian and a Poissonian distribution if both 
barriers have equal height and reduces to a Poissonian if one 
barrier is much higher than the other. 

PACS numbers: 72.70.+m, 72.10.Bg, 73.23.Ps, 02.50.Ey 



I. INTRODUCTION 

The nature of the current flow at low temperatures 
through mesoscopic structures has received a lot of at- 
tention during the last years. After initial focus on the 
conductance, which measures the average number of elec- 
trons transmitted in time, there has been an increasing 
interest in the noise power, a measure for the variance 
of the transmitted charge. At zero temperature, these 
current fluctuations are due to the discreteness of the 
electron charge. It has been found theoretically, that 
the zero-frequency shot-noise power P can be suppressed 
below its classical value characteristic for uncorrelated 
electron transport, -Ppoisson = 2el, with / the average 
current. This suppression is due to correlated electron 
transmission imposed by the Pauli principleEJtj. Conse- 
quently, it has been shown that in a double=.baxrier junc- 
tion P can be suppressed down to iPp i sson l3~Ej, depend- 
ing on the relative height of the barriers. For a metal- 
lic, diffusive eand-uctor various calculations yield that 
P = iPp oisson ll3ilZl. The shot-noise suppressio.D_.in_ these 
two systems has been observed experimentallyll31__. 

Recently, Levitov and Lesovik have gone one step fur- 
ther by studying the full distribution function roi charge 
transmitted through a mesoscopic conductor!--. This 
function gives the probability that a certain number of 
electrons are transmitted during a given time interval. 
Their quantum mechanical analysis demonstrates, that 
the attempts to transmit electrons are periodic in time, 
yielding a binomial distribution of transmitted electrons. 
On the basis of this result, Lee, Levitov, and Yakovets 
have calculated the charge distribution function for trans- 
port through a metallic, diffusive conductors. 



In this paper, we derive the complete distribution of 
transmitted charge through a double-barrier junction, by 
two different methods: Firstly, we follow a semiclassical 
approach, in which phase is neglected but the Pauli prin- 
ciple is accounted for. Here, the electron transport is 
described by classical jump rates. Secondly, we take a 
quantum mechanical approach, where we average the re- 
sult of Ref. ^3| over the distribution of transmission prob- 
abilities through the double-barrier system. For com- 
putational convenience, our analysis is carried out for a 
single-channel conductor in the absence of spin degen- 
eracy. However, the results can directly be generalized 
to a multi-channel conductor. Furthermore, we neglect 
charging effects, we assume small applied voltage as well 
as zero temperature, and we restrict ourselves to high 
tunnel barriers. 

The linear-response conductance G of a double-barrier 
junction is given by 



G = 



rir 2 



h Ti 



(1) 



with Ti <C 1 the transmission probability through barrier 
i = 1,2. Interestingly, the two approaches-tp derive Eq. 
(|l|) are of a completely different natureOEl The semi- 
classical derivation consists essentially in the addition of 
the resistances of both junctions, whereas the quantum 
mechanical derivation involves an average over resonant 
and non-resonant states. Physically, this averaging may 
correspond either to an applied voltage larger than the 
width of the resonance or to a summation over the modes 
in a multi-channel conductor if the distance between the 
barriers is larger than the Fermi wave length. 

The role of the presence of phase coherence on thcJjluep 
tuations in the current is still an intriguing issueOca. 
For example, the one-third suppression of the shot noise 
in a metallic, diffusive conductor was [Originally surmised 
to be of quantum mechanical origin__. However, later 
derivations through a semiclassical_aporoach yielded a 
suppression by one- third as welOEM With respect 
to the shot-noise power in the double-barrier junction, a 
quantum mechanical theory by Chen, and Tingcl and a 
semiclassical theory by Davies et alu give identical re- 
sults, namely 



P = 



iV + r2 

(r 1 + r 2 ; 



r^Poisson 



(2) 



An additional aim of the present paper is to check to 
which extent this insensitivity to the presence of phase 
coherence applies also for the complete distribution of 
transmitted charge. 
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The quantity of interest is P n (t), denoting the proba- 
bility that exactly n electrons have been transmitted dur- 
ing a time interval t. An alternative way to describe this 
distribution is through its characteristic-iunction x(A,i). 
They are mutually related according tcEa 



7T 

Pn(t) = ±- I d\e~ m \(\,t) 



(3a) 
(3b) 



It is often more convenient to determine x(A, t) instead of 
P n (t). Normalization requires that x(0, £) = 1< Further- 
more, it follows from Eq. (3a) that the average number 
of electrons transmitted during a time t is given by 







i(t) = nP n (t) = lim — X (A, t) 



(4) 



n=0 



More generally, one can express the fcth moment (t) of 
the distribution according to 



Hk(t) = n k (t) = lim 



d 



x{\t) 



(5) 



We note that the average current is simply / = efj,\(t)/t, 
whereas the noise power is proportional to the vari- 
ance of the number of transmitted electrons P = 
2e 2 lim t _ 00 varra(t)/* = 2e 2 lim t ^ 00 [ / u 2 (t) - m 2 (t)}/t. 
The logarithm of the characteristic function can be ex- 
panded as follows 



lnX(A,i)=2^/s fc (M 



fc=i 



(6) 



with Kk(t) the fcth cumulant of the distribution. The 
moment*, and the cumulants have a direct polynomial 
relationE3, for example fti(i) — Hi(t), K2(t) = ^{t) — 
Ml 2 (t). 

The .quantum mechanical analysis by Levitov and 
LesovikO has yielded the characteristic function of the 
charge through a single-channel conductor at zero tem- 
perature and at small voltage V , 



X(A,i) = [(e 



iX 



l)T 



leVt/h 



(7) 



with T the transmission probability at the Fermi level 
through the conductor. From Eqs. (||) and (Q) one can 
immediately derive the Landauer formula for the canduc- 
tance G and the formula for the shot-noise powemcl 



G = T T ' 



P 



2e 3 V 



T(l-T) 



(8) 
(9) 



Note that Eq. (g) is only valid in the phase-coherent 
regime, whereas in the absence of phase coherence P is 



given by a different equation!^ In Sec. |T| it is demon- 
strated how the charge distribution through a single bar- 
rier can be derived in a semiclassical approach. Sec. Ill 
repeats this analysis for the double-barrier junction. The 
quantum mechanical calculation is given in Sec. |jv|, after 
which we conclude in Sec. [v| 



II. SIMPLE EXAMPLE: SINGLE-BARRIER 
JUNCTION 

Let us illustrate our approach, by calculating the dis- 
tribution of transmitted charge through a single-channel, 
single-barrier junction, with a transmission probability 
r < 1. The average current through the barrier I = 
e 2 VT/h = e 7 , with 7 = eVY/h the tunnel rate through 
the barrier. The probability P n {t) that n electrons have 
been transmitted in a time t obeys the master equation 



dP (t) 
dt 

dP n {t) 

dt 



= -7P0W , (10a) 
= 7P n _i(t)- 7 P n (t), ifn>l, (10b) 



with the initial condition P„(0) = 5o,n- Eq. ( [io| ) can 
be solved straightforwardly by various means. Here, we 
adopt an approach, which appears to be useful for the 
double-barrier junction. The solution of Eq. (10a) is 
Po(t) = exp(— 7 i). We define the waiting-time distribu- 
tion t(>(t) = —dPo(t)/dt, denoting the probability density 
that an electron is transmitted immediately after having 
waited a time t, 



We write 



P n (t) = G n (t) - G n+1 (t) , 



(11) 



(12) 



where G n (t) denotes the probability that n or more elec- 
trons are transmitted during a time t. It can be calcu- 
lated according to 

t t t 
G n {t) = j dh J "dt a • • ■ J dt n 

tl tn-l 

x ip{ti)i){t 2 -ti)--' iP(t n - t n -x) . (13) 

Since this is a convolution, one has for the Laplace trans- 
form 



G n {s) = J dte- st G n (t) = - s l4,( s )] n , (14) 


where the Laplace transform of Eq. (nTl) is given by 



7 



s + 7 



(15) 
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From Eqs. (£|), (|T|), and (|T|), we find 



Pn(s) 



7" 



(s + 7) n+1 ' 
yielding for the distribution in time 



Pn(t) = 



-7* 



(16) 



(17) 



This is the Poisson distribution, as one would expect for 
the uncorrelated electron transfers through the barrier. 
The characteristic function is given by 



x(A,£)=exp[ 7 t(e* A -l)] 



(18) 



Indeed, in the limit T = Y -> 0, Eqs. (@) and coin- 
cide. 



III. CLASSICAL APPROACH 

We now study the double-barrier junction. The tunnel 
rate through barrier i = 1,2 is 7j = eVTi/h. Due to the 
Pauli principle, the number of electrons in the double- 
barrier junction can be either or 1. 

Our analysis is similar to the single-barrier case. How- 
ever, one now has to take into account two possible initial 
conditions at t = 0: either electrons in the junction — a 
situation with probability 72/(71 +72) — or 1 electron — 
probability 71/(71 +72)- The distribution of transmitted 
charge is thus given by 

P n (t) = _J_pW (t ) + _L P W (i ) , (19) 
71 +72 71 + 72 

where P„ (t) starts from j electrons in the junction at 
t = 0. The probability that at least n electrons have been 
transmitted can be expressed as 

G^(s) = -[Ms)Ms)] n , (20a) 
s 

G£Xs) = ±Ms)[M8)M8)] n - 1 , (20b) 

with &(«) = 7i/0 + 7*)- From Eqs. ©, ©, and © 
we obtain 

5 riA 5 (7i + 72) + 7i 2 + 72 2 + 7i72 , Q1 s 

^o(s) = — -. ■ w . T7— \ , (21a) 

(71 +72)(s + 7i)(s + 7 2 ) 

(7i72) n (g + 7i +72) 2 
(71 +72)(s + 7i) ,l+1 (s + 72) n+1 ' 
ifn>l. (21b) 



Pn(s) = 



The distribution function in time can be obtained 
through the inverse Laplace transform. In general, it 
yields a rather cumbersome expression. However, for the 
case of a symmetric double-barrier junction, 71 = 72 = 7, 
one has 



Po(*)=(l + |)e-^, (22a) 



Pn(t) 



2(2ri-l)! (2n)! 2(2n+l)! 
if n > 1 . 



(22b) 



For arbitrary 71 and 72, we can evaluate from Eq. ( |2l| ) 
the Laplace transform of the characteristic function 



1 



7i + 72 



(s + 71 + 7 2 ) 2 



(s + 7i)(s + 72) - e iA 7i 72 



1 



(23) 



yielding for the characteristic function in time 



X (X,t) = exp[-i( 7l + 72 )t] |cosh[|/3(A)t] 

27i72(^-in sinh[ ^ , (V)/ 



/3(A) 



7i+72 (7i+72)/?(A) 



(24) 



with /3(A) = v/( 7 i + 72 ) 2 +47 1 7 2 (e lA - 1). Eq. (g) is 
the central result of our semiclassical analysis. Let us 
evaluate from Eqs. (||) and (24) the first two moments of 
the charge distribution. The average number of electrons 
transmitted during a time t is 



n(t) = m(t) = 7l , 72 t , 
71 + 72 



(25) 



in agreement with Eq. (|]). For the second moment we 
find 

277\ M (7172) 2 2 , 7i72(7i 2 + 72 2 ) ^ 
n z (t) = H2(t) = 7 r^t H — t 

(71 +72) 2 (71 +72) 3 



2(7i7 2 ) 2 
(71 + 72) 4 



{l-exp[-(7i+ 72 )t]}. (26) 



The Fano factor, defined as the ratio of the variance 
to the average number of transmitted electrons, r(t) = 
[M*) - ^i 2 (i)]/^i(0> follows from Eqs. (pi) and <M), 



"(*) = 



7i 2 + 72 2 
(71 +72) 2 



27i72 (1 



-(71+72)^ 



(71 + J2) 3 t 



(27) 



The Fano factor gives the relative magnitude of the cur- 
rent fluctuations. Indeed, for large t, r(t) yields the shot- 
noise suppression according to Eqs. (Q) and (||). In Fig. |l| 
we have plotted r(t), for 71 = 72 = 7. We find that r(t) 
goes from 1 at small t, indicative for uncorrelated elec- 
tron transmission, to i at large t, indicative for a more 
correlated electron transmission. It is already within one 
percent of its final value at jt = 50, which corresponds 
to an average number of 25 transmitted electrons. 

The cumulants of the charge distribution can be deter- 
mined from the logarithm of Eq. (|2~i|) 
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ln X (V) = i[/3(A)-( 7 i+72)]t 

+ ln(l + e-^*) - In 2 

g(A) 2 7 i72(e' A - 1) 
7i+72 (7i+72)/?(A) 



In 



x tanh[i/3(A)t] 



(28) 



In general, it is cumbersome to derive the cumulants from 
this result. However, for large t, only the first term re- 
mains of importance. 



Substituting this result into Eqs. (||) and we recover 
for (G) and (P) the expressions given by Eqs. ([!]) and 

(D- 

In order to obtain the ensemble average of all the cu- 
mulants of the distribution function, _ss£e average the log- 
arithm of the characteristic functioned. For the double- 
barrier junction, we obtain from Eq. (ffl) 



eVt 

(lnx(A,i)) = Po— I dv In 



(e lX - l)?i 
1 + v 2 



(32) 



IV. QUANTUM MECHANICAL APPROACH 

Whereas in a semiclassical picture, the transmission 
probability through the double-barrier junction is just a 
constant, in a quantum mechanical approach, the trans- 
mission probability T varies according to a Fabry-Perot 
type of formula 



T 



r,r. 



11 2 



i - 2^(1 - fT)(i - r 2 ) cos0 + (i - ri)(i - r 2 ) ' 

(29) 

where <fi is the phase accumulated in one round trip be- 
tween the barriers. The distribution function p(T) of 
the transmission probabilities through the system can be 
obtained from the assumption that <j> is uniformly dis- 
tributed between and 2ir£3 In the limit I\, T 2 <C 1 this 
implies 



P{T) 



rir 2 



i 



7r(Ti +r 2 ) ^/T 3 (Ti — T) 



(30) 



if T G [T ,Ti], and p(T) = otherwise, with T = 

4r 1 r 2 /[(r 1 + r 2 ) 2 + 4^ 2 ] and t x = 4r 1 r 2 /(r 1 + r 2 ) 2 . 

The distribution function (|3(]) is plotted in Fig. ||. Simi- 
lar to a metallic, diffusive conductort 3 !, this distribution 
is bimodal, in the sense that the transmission probabili- 
ties are either close to To ~ or close to T\. 

The ensemble average of a quantity a(T) over all 
possible transmission probabilities is given by (a) = 

J^ 1 dTp(T)a(T). It is convenient to switch variables from 
T to v with T = Ti/(1 + v 2 ), so that p(v) = p = 
(Tx +r 2 )/27r is uniform over the range [0, u max \. In prac- 
tice, the upper limit v max can be often replaced by infin- 
ity. The ensemble average for the mth power (m > 1) of 
the transmission probability is given by 



(T m ) 



1 



Po I 'I'' 

o 

(2m -2)! 

[(m-i)!] 2 (r 1 + r 2 )^ 



(rir 2 



(31) 



with /3(A) = \/(7i + 72) 2 + 47i7 2 (e lA - 1). This is the 
key result of the quantum mechanical evaluation. Using 
Eq. (||) , we find for the ensemble average of the first three 
cumulants 



<«!(*)> 



7i 72 
7i + 72 



t . 



, f.s, 7i72(7i 2 +72 2 ) , 
{^2{t)) = — - — — — — * , 



(«s(*)> = 



(7i +72) 3 
7i 72 



(33a) 
(33b) 



(71 + 72) 5 



x ( 7l 4 - 2 7 i 3 72 + 67i 2 7 2 2 - 2 7 i7 2 3 + 72 4 ) . 

(33c) 

Since (ln%(A,t)) is proportional to t, all the cumulants 
are linear in t as well. This implies for the Fano factor 



(r(t)) = {K2{t)) 



71 



72 



(«i(0> (71 + 72) 2 



(34) 



which is constant in time and equal to the large-i value of 
Eq. @ (see Fig. |). If 71 > 72 (or vice versa), Eq. @ 
reduces to a Poissonian distribution, as expected. For a 



symmetric double-barrier junction with 7! 
expression (|32| ) simplifies considerably: 

(ln X (A,<))= 7 i(e lA/2 -l) . 

For the feth cumulant we find from Eqs. (^ 

7 t 



7 2 = 7, the 
(35) 



and (B5 



(«fc(*)> = 



2 k ' 



(36) 



The charge distribution is thus somewhere between Gaus- 
sian (where Kk — for k > 3) and Poissonian (where 
Kk = it for all k) . 



V. CONCLUSIONS 

Let us make the comparison between the outcome of 
the two approaches, i.e., between the semiclassical re- 
sult ln%(A,i) from Eq. (|28| ) and the quantum result 
(m%(A, t)) given in Eq. (|32j), where we have averaged 
over the distribution of transmission probabilities. The 
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semiclassical result yields a more complicated expression, 
however the most important contribution, which is pro- 
portional to t, is precisely equivalent to the result (|32|). 
The other terms are either constant (do not depend on 
t) or vanish exponentially with t. Only on short time 
scales, corresponding to the transfer of a few electrons, 
we see sizeable differences between both approaches. We 
surmise that these differences at small t are not due to 
the neglect of the phase, but rather depend on the pre- 
cise way the reservoirs are modeled in both approaches. 
The quantum mechanical derivation is based on Eq. (7j), 
which applies for arbitrary transmission probability be- 
tween and 1. Here, the number of electrons transmitted 
has a maximum value n(t) = eVt/h. Our semiclassi- 
cal calculation assumes independent tunnel events and is 
therefore only valid for small transmission probabilities. 
However, the number of electrons which can be transmit- 
ted n(t) is not bounded. Even though this seems not to 
be very important for the case of high tunnel barriers, 
one may expect that it leads to differences on a small 
time scale. Therefore, we just draw conclusions from the 
comparison at larger t. Here we find that the two results 
are equal, and that as a consequence, all the cumulants 
are also equal. This demonstrates that the statistics of 
charge transport through a double-barrier junction does 
not reveal whether phase coherence is present or absent. 

This inscnsitivity to the presence of phase coherence, 
does not imply that phase breaking is not of influence 
in a real experiment. This depends on the physical pro- 
cess which destroys the phase coherence. For example, 
using the method given in Ref. [HJ one can show that 
electron-electron scattering inside the double-barrier sys- 
tem, in which both phase coherence is destroyed and en- 
ergy is redistributed among the electrons, increases the 
shot noise above the value of Eq. (^). However, the analy- 
sis in the present paper demonstrates that merely break- 
ing the phase leaves the charge transport through the 
system unaffected. This is in contrast to the result of 
Ref. in which incoherence is modeled by adding ran- 
dom phases to the wave function on each round trip. The 
authors find that this increases the shot noise, so that we 
conclude that their model is not equivalent to just de- 
stroying the phase. 

It might be interesting to determine the role of inelastic 
processes inside the tunnel barrier on the charge distri- 
bution. We think that these effects can well be taken into 
account using the semiclassical analysis, whereas a com- 
plete quantum mechanical derivation looks more compli- 
cated. Another extension of the work described in this 
paper, would be to repeat the semiclassical analysis for a 
metallic, diffusive conductor, and compare the outcome 
with the quantum mechanical derivation of Ref. |24|. 

In summary, we have derived the complete distribution 
of transmitted charge through a double-barrier junction 
at zero temperature and at low voltage. We have used 
a semiclassical approach on the basis of classical jump 
rates as well as a quantum mechanical approach, in which 
the result of Levitov and Lesovik for an arbitrary single- 



channel conductoill is averaged over the distribution of 
transmission probabilities through the system. Our re- 
sults are in precise agreement with previous, values for 
the conductance and for the shot-noise powerQ'Q. Within 
both approaches, we have determined the logarithm of 
the characteristic function, which become equivalent at 
large times. It is found that for symmetric tunnel barri- 
ers, the charge distribution is between a Gaussian and a 
Poissonian distribution. 
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FIG. 2. The distribution of transmission probabilities 
through a double-barrier junction in a quantum mechanical 
model, according to Eq. @, for Ti = 0.02 and T 2 = 0.03. 
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FIG. 1. The Fano factor r(t), giving the ratio of the vari- 
ance to the average number of transmitted electrons, versus 
time t for a symmetric double-barrier junction with tunnel 
rates 71 = 72 = 7. The solid line gives the result (^) of the 
semiclassical analysis and the dashed line the quantum result 
(r(i)), according to Eq. (fjij). 
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